Genome-wide association analysis of composite sleep health scores in 413,904 individuals.

Recent genome-wide association studies (GWASs) of several individual sleep traits have identified hundreds of genetic loci, suggesting diverse mechanisms. Moreover, sleep traits are moderately correlated, and together may provide a more complete picture of sleep health, while also illuminating distinct domains. Here we construct novel sleep health scores (SHSs) incorporating five core self-report measures: sleep duration, insomnia symptoms, chronotype, snoring, and daytime sleepiness, using additive (SHS-ADD) and five principal components-based (SHS-PCs) approaches. GWASs of these six SHSs identify 28 significant novel loci adjusting for multiple testing on six traits (p<8.3e-9), along with 341 previously reported loci (p<5e-08). The heritability of the first three SHS-PCs equals or exceeds that of SHS-ADD (SNP-h2=0.094), while revealing sleep-domain-specific genetic discoveries. Significant loci enrich in multiple brain tissues and in metabolic and neuronal pathways. Post GWAS analyses uncover novel genetic mechanisms underlying sleep health and reveal connections to behavioral, psychological, and cardiometabolic traits.

Heritability and genetic correlations for sleep apnea, insomnia, and hypersomnia in a large clinical biobank.

RATIONALE: Comorbid insomnia and sleep apnea is reported to have worse outcomes than either condition alone. The local genetic correlations of these disorders are unknown. OBJECTIVES: To identify local genomic regions with heritability for clinically diagnosed sleep apnea and insomnia, and to identify local genetic correlations between these disorders and/or hypersomnia. METHODS: Fifty thousand two hundred seventeen patients of European ancestry were examined. Global and local heritability and genetic correlations for independent regions were calculated, adjusting for obesity and other covariates. RESULTS: Sleep apnea and insomnia were significantly globally heritable and had 118 and 168 genetic regions with local heritability p-values <.05, respectively. One region had a significant genetic correlation for sleep apnea and hypersomnia (p-value = 9.85 × 10-4). CONCLUSIONS: Clinically diagnosed sleep apnea and insomnia have minimal shared genetic architecture, supporting genetically distinct comorbid insomnia and sleep apnea components. However, additional correlated regions may be identified with additional sample size and methodological improvements.

A statistical framework for powerful multi-trait rare variant analysis in large-scale whole-genome sequencing studies.

Large-scale whole-genome sequencing (WGS) studies have improved our understanding of the contributions of coding and noncoding rare variants to complex human traits. Leveraging association effect sizes across multiple traits in WGS rare variant association analysis can improve statistical power over single-trait analysis, and also detect pleiotropic genes and regions. Existing multi-trait methods have limited ability to perform rare variant analysis of large-scale WGS data. We propose MultiSTAAR, a statistical framework and computationally-scalable analytical pipeline for functionally-informed multi-trait rare variant analysis in large-scale WGS studies. MultiSTAAR accounts for relatedness, population structure and correlation among phenotypes by jointly analyzing multiple traits, and further empowers rare variant association analysis by incorporating multiple functional annotations. We applied MultiSTAAR to jointly analyze three lipid traits (low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides) in 61,861 multi-ethnic samples from the Trans-Omics for Precision Medicine (TOPMed) Program. We discovered new associations with lipid traits missed by single-trait analysis, including rare variants within an enhancer of NIPSNAP3A and an intergenic region on chromosome 1.

Characterization of sleep apnea among a sample of adults from Samoa.

Sleep apnea is a public health concern around the world, but little research has been dedicated to examining this issue in low- and middle-income countries, including Samoa. Using data collected through the Soifua Manuia ("Good Health") study, which aimed to investigate the impact of the body mass index (BMI)-associated genetic variant rs373863828 in CREB3 Regulatory Factor ( CREBRF ) on metabolic traits in Samoan adults, we examined the sample prevalence and characteristics of sleep apnea using data collected with a validated home sleep apnea device (WatchPAT, Itamar). A total of 330 participants (sampled to overrepresent the obesity-risk allele of interest) had sleep data available. Participants (53.3% female) had a mean (SD) age of 52.0 (9.9) years and BMI of 35.5 (7.5) kg/m 2 and 36.3% of the sample had type 2 diabetes. Based on the 3% and 4% apnea hypopnea indices (AHI) and the 4% oxygen desaturation index (ODI), descriptive analyses revealed that many participants had potentially actionable sleep apnea defined as >5 events/hr (87.9%, 68.5%, and 71.2%, respectively) or clinically actionable sleep apnea defined as ≥15 events/hr (54.9%, 31.5%, and 34.5%, respectively). Sleep apnea was more severe in men; for example, clinically actionable sleep apnea (≥15) based on the AHI 3% definition was observed in 61.7% of men and 48.9% of women. Correction for non-representational sampling related to the CREBRF obesity-risk allele resulted in only slightly lower estimates. Across the AHI 3%, AHI 4%, and ODI 4%, multiple linear regression revealed associations between a greater number of events/hr and higher age, male sex, higher body mass index, higher abdominal-hip circumference ratio, and geographic region of residence. Our study identified a much higher frequency of sleep apnea in Samoa compared with published data from other studies, but similar predictors. Continued research addressing generalizability of these findings, as well as a specific focus on diagnosis and affordable and equitable access to treatment, is needed to alleviate the burden of sleep apnea in Samoa and around the world.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed whole-genome sequencing study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions in lipid metabolism. Large-scale whole-genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess more associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with measurement of blood lipids and lipoproteins (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare-variant aggregate association tests using the STAAR (variant-set test for association using annotation information) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare-coding variants in nearby protein-coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500-kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variation and rare protein-coding variation at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNAs.

Rare variants in long non-coding RNAs are associated with blood lipid levels in the TOPMed Whole Genome Sequencing Study.

Long non-coding RNAs (lncRNAs) are known to perform important regulatory functions. Large-scale whole genome sequencing (WGS) studies and new statistical methods for variant set tests now provide an opportunity to assess the associations between rare variants in lncRNA genes and complex traits across the genome. In this study, we used high-coverage WGS from 66,329 participants of diverse ancestries with blood lipid levels (LDL-C, HDL-C, TC, and TG) in the National Heart, Lung, and Blood Institute (NHLBI) Trans-Omics for Precision Medicine (TOPMed) program to investigate the role of lncRNAs in lipid variability. We aggregated rare variants for 165,375 lncRNA genes based on their genomic locations and conducted rare variant aggregate association tests using the STAAR (variant-Set Test for Association using Annotation infoRmation) framework. We performed STAAR conditional analysis adjusting for common variants in known lipid GWAS loci and rare coding variants in nearby protein coding genes. Our analyses revealed 83 rare lncRNA variant sets significantly associated with blood lipid levels, all of which were located in known lipid GWAS loci (in a ±500 kb window of a Global Lipids Genetics Consortium index variant). Notably, 61 out of 83 signals (73%) were conditionally independent of common regulatory variations and rare protein coding variations at the same loci. We replicated 34 out of 61 (56%) conditionally independent associations using the independent UK Biobank WGS data. Our results expand the genetic architecture of blood lipids to rare variants in lncRNA, implicating new therapeutic opportunities.

The public health impact of poor sleep on severe COVID-19, influenza and upper respiratory infections.

BACKGROUND: Poor sleep is associated with an increased risk of infections and all-cause mortality but the causal direction between poor sleep and respiratory infections has remained unclear. We examined if poor sleep contributes as a causal risk factor to respiratory infections. METHODS: We used data on insomnia, influenza and upper respiratory infections (URIs) from primary care and hospital records in the UK Biobank (N ≈ 231,000) and FinnGen (N ≈ 392,000). We computed logistic regression to assess association between poor sleep and infections, disease free survival hazard ratios, and performed Mendelian randomization analyses to assess causality. FINDINGS: Utilizing 23 years of registry data and follow-up, we discovered that insomnia diagnosis associated with increased risk for infections (FinnGen influenza Cox's proportional hazard (CPH) HR = 4.34 [3.90, 4.83], P = 4.16 × 10-159, UK Biobank influenza CPH HR = 1.54 [1.37, 1.73], P = 2.49 × 10-13). Mendelian randomization indicated that insomnia causally predisposed to influenza (inverse-variance weighted (IVW) OR = 1.65, P = 5.86 × 10-7), URI (IVW OR = 1.94, P = 8.14 × 10-31), COVID-19 infection (IVW OR = 1.08, P = 0.037) and risk of hospitalization from COVID-19 (IVW OR = 1.47, P = 4.96 × 10-5). INTERPRETATION: Our findings indicate that chronic poor sleep is a causal risk factor for contracting respiratory infections, and in addition contributes to the severity of respiratory infections. These findings highlight the role of sleep in maintaining sufficient immune response against pathogens. FUNDING: Instrumentarium Science Foundation, Academy of Finland, Signe and Ane Gyllenberg Foundation, National Institutes of Health.

The genetic determinants of recurrent somatic mutations in 43,693 blood genomes.

Nononcogenic somatic mutations are thought to be uncommon and inconsequential. To test this, we analyzed 43,693 National Heart, Lung and Blood Institute Trans-Omics for Precision Medicine blood whole genomes from 37 cohorts and identified 7131 non-missense somatic mutations that are recurrently mutated in at least 50 individuals. These recurrent non-missense somatic mutations (RNMSMs) are not clearly explained by other clonal phenomena such as clonal hematopoiesis. RNMSM prevalence increased with age, with an average 50-year-old having 27 RNMSMs. Inherited germline variation associated with RNMSM acquisition. These variants were found in genes involved in adaptive immune function, proinflammatory cytokine production, and lymphoid lineage commitment. In addition, the presence of eight specific RNMSMs associated with blood cell traits at effect sizes comparable to Mendelian genetic mutations. Overall, we found that somatic mutations in blood are an unexpectedly common phenomenon with ancestry-specific determinants and human health consequences.

Aberrant activation of TCL1A promotes stem cell expansion in clonal haematopoiesis.

Mutations in a diverse set of driver genes increase the fitness of haematopoietic stem cells (HSCs), leading to clonal haematopoiesis1. These lesions are precursors for blood cancers2-6, but the basis of their fitness advantage remains largely unknown, partly owing to a paucity of large cohorts in which the clonal expansion rate has been assessed by longitudinal sampling. Here, to circumvent this limitation, we developed a method to infer the expansion rate from data from a single time point. We applied this method to 5,071 people with clonal haematopoiesis. A genome-wide association study revealed that a common inherited polymorphism in the TCL1A promoter was associated with a slower expansion rate in clonal haematopoiesis overall, but the effect varied by driver gene. Those carrying this protective allele exhibited markedly reduced growth rates or prevalence of clones with driver mutations in TET2, ASXL1, SF3B1 and SRSF2, but this effect was not seen in clones with driver mutations in DNMT3A. TCL1A was not expressed in normal or DNMT3A-mutated HSCs, but the introduction of mutations in TET2 or ASXL1 led to the expression of TCL1A protein and the expansion of HSCs in vitro. The protective allele restricted TCL1A expression and expansion of mutant HSCs, as did experimental knockdown of TCL1A expression. Forced expression of TCL1A promoted the expansion of human HSCs in vitro and mouse HSCs in vivo. Our results indicate that the fitness advantage of several commonly mutated driver genes in clonal haematopoiesis may be mediated by TCL1A activation.

Impact of Pulmonary and Sleep Disorders on COVID-19 Infection Severity in a Large Clinical Biobank.

Rationale: Multiple pulmonary, sleep, and other disorders are associated with the severity of Covid-19 infections but may or may not directly affect the etiology of acute Covid-19 infection. Identifying the relative importance of concurrent risk factors may prioritize respiratory disease outbreaks research. Objectives: To identify associations of common preexisting pulmonary and sleep disease on acute Covid-19 infection severity, investigate the relative contributions of each disease and selected risk factors, identify sex-specific effects, and examine whether additional electronic health record (EHR) information would affect these associations. Methods: 45 pulmonary and 6 sleep diseases were examined in 37,020 patients with Covid-19. We analyzed three outcomes: death; a composite measure of mechanical ventilation and/or ICU admission; and inpatient admission. The relative contribution of pre-infection covariates including other diseases, laboratory tests, clinical procedures, and clinical note terms was calculated using LASSO. Each pulmonary/sleep disease model was then further adjusted for covariates. Measurements and main results: 37 pulmonary/sleep diseases were associated with at least one outcome at Bonferroni significance, 6 of which had increased relative risk in LASSO analyses. Multiple prospectively collected non-pulmonary/sleep diseases, EHR terms and laboratory results attenuated the associations between preexisting disease and Covid-19 infection severity. Adjustment for counts of prior "blood urea nitrogen" phrases in clinical notes attenuated the odds ratio point estimates of 12 pulmonary disease associations with death in women by ≥1. Conclusions: Pulmonary diseases are commonly associated with Covid-19 infection severity. Associations are partially attenuated by prospectively-collected EHR data, which may aid in risk stratification and physiological studies.

Obstructive Sleep Apnea and Longitudinal Changes in Interstitial Lung Imaging and Lung Function: The MESA Study.

Rationale: Obstructive sleep apnea (OSA) has been hypothesized to be a risk factor in interstitial lung disease (ILD) and is associated with radiological markers that may represent the earlier stages of ILD. Prior studies have been limited by their cross-sectional design and potential confounding by body habitus. Objectives: To test the hypothesis that OSA severity is associated with more high-attenuation areas (HAAs) on computed tomography and worse lung function over time among older community-dwelling adults. Methods: We used data from participants in the MESA (Multi-Ethnic Study of Atherosclerosis) who had apnea-hypopnea index (AHI) measured from polysomnography (2010-2013), high attenuation areas (HAAs, -600 to -250 Hounsfield units, n = 784), assessments from exams 5 (2010-2012) and 6 (2016-2018) full-lung computed tomography scans, and spirometry assessments (n = 677). Linear mixed-effects models with random intercept were used to examine associations of OSA severity (i.e., AHI and hypoxic burden) with changes in HAAs, total lung volumes, and forced vital capacity (FVC) between exams 5 and 6. Potential confounders were adjusted for in the model, including age, sex, smoking history, height, and weight. Results: Among those with a higher AHI there were more men and a higher body mass index. Participants with AHI ⩾ 15 events/h and in the highest hypoxic burden quartile each had increases in HAAs of 11.30% (95% confidence interval [CI], 3.74-19.35%) and 9.85% (95% CI, 1.40-19.01%) per 10 years, respectively. There was a more rapid decline in total lung volumes imaged and FVC among those with AHI ⩾ 15 events/h of 220.2 ml (95% CI, 47.8-392.5 ml) and 3.63% (95% CI, 0.43-6.83%) per 10 years, respectively. Conclusions: A greater burden of hypoxia related to obstructive events during sleep was associated with increased lung densities over time and a more rapid decline in lung volumes regardless of body habitus. Our findings suggest OSA may be a contributing factor in the early stages of ILD.

Structural variation across 138,134 samples in the TOPMed consortium.

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hematologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

Structural variation across 138,134 samples in the TOPMed consortium.

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hemotologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

Powerful, scalable and resource-efficient meta-analysis of rare variant associations in large whole genome sequencing studies.

Meta-analysis of whole genome sequencing/whole exome sequencing (WGS/WES) studies provides an attractive solution to the problem of collecting large sample sizes for discovering rare variants associated with complex phenotypes. Existing rare variant meta-analysis approaches are not scalable to biobank-scale WGS data. Here we present MetaSTAAR, a powerful and resource-efficient rare variant meta-analysis framework for large-scale WGS/WES studies. MetaSTAAR accounts for relatedness and population structure, can analyze both quantitative and dichotomous traits and boosts the power of rare variant tests by incorporating multiple variant functional annotations. Through meta-analysis of four lipid traits in 30,138 ancestrally diverse samples from 14 studies of the Trans Omics for Precision Medicine (TOPMed) Program, we show that MetaSTAAR performs rare variant meta-analysis at scale and produces results comparable to using pooled data. Additionally, we identified several conditionally significant rare variant associations with lipid traits. We further demonstrate that MetaSTAAR is scalable to biobank-scale cohorts through meta-analysis of TOPMed WGS data and UK Biobank WES data of ~200,000 samples.

Whole genome sequence analysis of apparent treatment resistant hypertension status in participants from the Trans-Omics for Precision Medicine program.

Introduction: Apparent treatment-resistant hypertension (aTRH) is characterized by the use of four or more antihypertensive (AHT) classes to achieve blood pressure (BP) control. In the current study, we conducted single-variant and gene-based analyses of aTRH among individuals from 12 Trans-Omics for Precision Medicine cohorts with whole-genome sequencing data. Methods: Cases were defined as individuals treated for hypertension (HTN) taking three different AHT classes, with average systolic BP ≥ 140 or diastolic BP ≥ 90 mmHg, or four or more medications regardless of BP (n = 1,705). A normotensive control group was defined as individuals with BP < 140/90 mmHg (n = 22,079), not on AHT medication. A second control group comprised individuals who were treatment responsive on one AHT medication with BP < 140/ 90 mmHg (n = 5,424). Logistic regression with kinship adjustment using the Scalable and Accurate Implementation of Generalized mixed models (SAIGE) was performed, adjusting for age, sex, and genetic ancestry. We assessed variants using SKAT-O in rare-variant analyses. Single-variant and gene-based tests were conducted in a pooled multi-ethnicity stratum, as well as self-reported ethnic/racial strata (European and African American). Results: One variant in the known HTN locus, KCNK3, was a top finding in the multi-ethnic analysis (p = 8.23E-07) for the normotensive control group [rs12476527, odds ratio (95% confidence interval) = 0.80 (0.74-0.88)]. This variant was replicated in the Vanderbilt University Medical Center's DNA repository data. Aggregate gene-based signals included the genes AGTPBP, MYL4, PDCD4, BBS9, ERG, and IER3. Discussion: Additional work validating these loci in larger, more diverse populations, is warranted to determine whether these regions influence the pathobiology of aTRH.

Leveraging pleiotropy to discover and interpret GWAS results for sleep-associated traits.

Genetic association studies of many heritable traits resulting from physiological testing often have modest sample sizes due to the cost and burden of the required phenotyping. This reduces statistical power and limits discovery of multiple genetic associations. We present a strategy to leverage pleiotropy between traits to both discover new loci and to provide mechanistic hypotheses of the underlying pathophysiology. Specifically, we combine a colocalization test with a locus-level test of pleiotropy. In simulations, we show that this approach is highly selective for identifying true pleiotropy driven by the same causative variant, thereby improves the chance to replicate the associations in underpowered validation cohorts and leads to higher interpretability. Here, as an exemplar, we use Obstructive Sleep Apnea (OSA), a common disorder diagnosed using overnight multi-channel physiological testing. We leverage pleiotropy with relevant cellular and cardio-metabolic phenotypes and gene expression traits to map new risk loci in an underpowered OSA GWAS. We identify several pleiotropic loci harboring suggestive associations to OSA and genome-wide significant associations to other traits, and show that their OSA association replicates in independent cohorts of diverse ancestries. By investigating pleiotropic loci, our strategy allows proposing new hypotheses about OSA pathobiology across many physiological layers. For example, we identify and replicate the pleiotropy across the plateletcrit, OSA and an eQTL of DNA primase subunit 1 (PRIM1) in immune cells. We find suggestive links between OSA, a measure of lung function (FEV1/FVC), and an eQTL of matrix metallopeptidase 15 (MMP15) in lung tissue. We also link a previously known genome-wide significant peak for OSA in the hexokinase 1 (HK1) locus to hematocrit and other red blood cell related traits. Thus, the analysis of pleiotropic associations has the potential to assemble diverse phenotypes into a chain of mechanistic hypotheses that provide insight into the pathogenesis of complex human diseases.

Genetic diversity fuels gene discovery for tobacco and alcohol use.

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury1-4. These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries5. Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction.

Correlates of daytime sleepiness and insomnia among adults in Samoa.

Objective: To describe daytime sleepiness and insomnia among adults in Samoa and identify modifiable factors associated with these measures. Design/setting: Cross-sectional analysis of data from the Soifua Manuia ("Good Health") study (n = 519, 55.1% female); Upolu island, Samoa. Methods: Daytime sleepiness and insomnia were assessed with the Epworth Sleepiness Scale (ESS) and the Women's Health Initiative Insomnia Rating Scale (WHIIRS), respectively. Detailed physical, sociodemographic, and behavioral factors were collected. Sleep measures were characterized using multiple linear regression with backwards elimination and a bootstrap stability investigation. Results: Excessive daytime sleepiness (ESS>10) and insomnia (WHIIRS>10) were reported by 20% and 6.3% of the sample, respectively. ESS scores were higher in individuals reporting more physical activity (Estimate=1.88; 95% CI=1.12 to 2.75), higher material wealth (0.18; 0.09 to 0.28), and asthma (2.85; 1.25 to 4.51). ESS scores were lower in individuals residing in periurban versus urban regions (-1.43; -2.39 to -0.41), reporting no work versus day shift work (-2.26; -3.07 to -1.41), and reporting greater perceived stress (-0.14; -0.23 to -0.06). WHIIRS scores were lower in individuals reporting "other" shift work (split/irregular/on-call/rotating) versus day shift work (-1.96; -2.89 to -1.14) and those who perceived their village's wealth to be poor/average versus wealthy (-0.94; -1.50 to -0.34). Conclusions: Participants had a generally higher prevalence of excessive daytime sleepiness, but lower prevalence of insomnia, compared with individuals from high-income countries. Factors associated with sleep health differed compared with prior studies, emphasizing potential cultural/environmental differences and the need for targeted interventions to improve sleep health in this setting.

A framework for detecting noncoding rare-variant associations of large-scale whole-genome sequencing studies.

Large-scale whole-genome sequencing studies have enabled analysis of noncoding rare-variant (RV) associations with complex human diseases and traits. Variant-set analysis is a powerful approach to study RV association. However, existing methods have limited ability in analyzing the noncoding genome. We propose a computationally efficient and robust noncoding RV association detection framework, STAARpipeline, to automatically annotate a whole-genome sequencing study and perform flexible noncoding RV association analysis, including gene-centric analysis and fixed window-based and dynamic window-based non-gene-centric analysis by incorporating variant functional annotations. In gene-centric analysis, STAARpipeline uses STAAR to group noncoding variants based on functional categories of genes and incorporate multiple functional annotations. In non-gene-centric analysis, STAARpipeline uses SCANG-STAAR to incorporate dynamic window sizes and multiple functional annotations. We apply STAARpipeline to identify noncoding RV sets associated with four lipid traits in 21,015 discovery samples from the Trans-Omics for Precision Medicine (TOPMed) program and replicate several of them in an additional 9,123 TOPMed samples. We also analyze five non-lipid TOPMed traits.

Whole genome sequence analysis of blood lipid levels in >66,000 individuals.

Blood lipids are heritable modifiable causal factors for coronary artery disease. Despite well-described monogenic and polygenic bases of dyslipidemia, limitations remain in discovery of lipid-associated alleles using whole genome sequencing (WGS), partly due to limited sample sizes, ancestral diversity, and interpretation of clinical significance. Among 66,329 ancestrally diverse (56% non-European) participants, we associate 428M variants from deep-coverage WGS with lipid levels; ~400M variants were not assessed in prior lipids genetic analyses. We find multiple lipid-related genes strongly associated with blood lipids through analysis of common and rare coding variants. We discover several associated rare non-coding variants, largely at Mendelian lipid genes. Notably, we observe rare LDLR intronic variants associated with markedly increased LDL-C, similar to rare LDLR exonic variants. In conclusion, we conducted a systematic whole genome scan for blood lipids expanding the alleles linked to lipids for multiple ancestries and characterize a clinically-relevant rare non-coding variant model for lipids.